Seeking Inflation Fossils in the Cosmic Microwave Background 
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If during inflation the inflaton couples to a "fossil" field, some new scalar, vector, or tensor field, 
it typically induces a scalar-scalar-fossil bispectrum. Even if the fossil field leaves no direct physical 
trace after inflation, it gives rise to correlations between different Fourier modes of the curvature 
or, equivalently, a nonzero curvature trispectrum. Here we quantify the effects of a fossil field on 
the cosmic microwave background (CMB) temperature fluctuations in terms of bipolar spherical 
harmonics (BiPoSHs). The effects of vector and tensor fossils can be distinguished geometrically 
from those of scalars through the parity of the BiPoSHs they induce. However, the two-dimensional 
nature of the CMB sky does not allow vectors to be distinguished geometrically from tensors. We 
estimate the detectability of a signal in terms of the scalar-scalar-fossil coupling for scalar, vector, 
and tensor fossils, assuming a local-type coupling. We comment on a divergence that arises in the 
quadrupolar BiPoSH from the scalar-scalar-tensor correlation in single-fleld slow-roll inflation. 



I. INTRODUCTION 

Despite the observational success of the inflation 
paradigm, the nature of the inflationary epoch is still 
poorly understood. The single-field slow-roll model, be- 
ing the simplest inflation model, might not be the whole 
story of the physics behind infiation. An abundance of 
inflation models introduce new fields, which are often 
coupled to the infiaton [l|, Q • Candidates for those new 
degrees of freedom include extra scalar fields in exten- 
sions of the simplest inflation model Q . They might also 
be primordial vector fields that drive or impact infla- 
tion Q , or even new vectorial degrees of freedom arising 
from modifled theories of gravity Even in the sim- 
plest single-fleld slow-roll inflation, tensor metric pertur- 
bations [y] couple to the inflaton jTHlOj. 

Fortunately, the cosmic microwave background (CMB) 
provides a special window to probe physics at very early 
moments, by encoding non-trivial primordial correlations 
beyond Gaussianity in its anisotropy pattern. By mea- 
suring the primordial bispectra, trispectra, and possibly 
higher-order correlation functions from its temperature 
and polarization anisotropics [lll-fisj. various inflation 
models will hopefully be distinguished, and new physics 
might be revealed. 

Recently, Ref. [l6j proposed a generic parameteriza- 
tion of the primordial bispectra involving the scalar met- 
ric perturbation $ (the gauge-invariant Bardeen poten- 
tial ) and a new field h from inflation. The new field is 
dubbed an inflation fossil, a hypothesized primordial de- 
gree of freedom that no longer interacts or very weakly 
interacts during late-time cosmic evolution. The only 
observational effect of an inflation fossil might therefore 
be its imprint in the primordial curvature perturbation. 
Inflation fossils can be those extra flelds that are intro- 
duced in a variety of alternatives to the single-fleld slow- 
roll model. Different models may be distinguished by the 
spin of the new fleld. To treat different possibilities for 
the spin in a model-independent fashion, in Ref. [l6j the 
fossil field is parameterized by a symmetric traceless ten- 



sor field. Since a symmetric traceless tensor contains a 
longitudinal scalar (L), two divergence- free vectors (V), 
and two divergence- free traceless tensors (T), the param- 
eterization allows for all three possibilities for the spin. 
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Figure 1: Primordial correlations of scalar perturbation (solid 
vectors) induced by the fossil field (dashed vectors) in Fourier 
space: (a) For a given realization of the fossil field, two differ- 
ent scalar-perturbation modes are correlated with each other, 
(b) For a stochastic background of fossil fields, a connected 
trispectrum of scalar perturbations is induced. 

A given realization of the fossil fleld causes two dif- 
ferent scalar perturbation modes to correlate with each 
other (panel (a) of Fig. [1]). Scalar, vector, and tensor fos- 
sils give rise to geometrically distinct cross-correlations. 
For a stochastic fossil background, connected trispectra 
that correlate four differents modes of $ are generated 
(panel (b) of Fig. [1]) when an ensemble average over all 
fossil-fleld realizations is performed. Such a four-point 
correlation can be sought in galaxy surveys where the 
correlations from scalar, vector and tensor fossils can be 
disentangled geometrically [l6j . 

In this paper we study the CMB signatures of fossil 
flelds. Correlations between different Fourier modes of 
the curvature perturbation induce couplings between dif- 
ferent CMB-temperature spherical-harmonic coefficients. 
We parametrize these cross-correlations in terms of bipo- 
lar spherical harmonics (BiPoSHs). In the presence of 
some realization of the fossil field, the CMB appears non- 
Gaussian with nonvanishing BiPoSHs, or equivalently, a 
nonzero four-point correlation function, or trispectrum. 
Here we first calculate the BiPoSH coefficients from a 



2 



specific realization of the fossil field, and then find the Bi- 
PoSH power spectra for a stochastic background for the 
fossil field. We write down the minimum-variance CMB 
estimators for the amplitudes of the fossil-field power 
spectra. As an example, we numerically examine the 
case of a local-type scalar-scalar-fossil bispectrum with a 
scale-free fossil spectrum. That case can be described 
by two parameters, the amplitude of the bispectrum, 
and the normalization of the fossil power spectrum. We 
find that the dominant effects from a local-type bispec- 
trum are modulations of small-angular-scale correlations 
on large angular scales. We evaluate the sensitivities for 
the reduced amplitude of the fossil background, a com- 
bination of the two parameters, for different fossil-field 
spin. Our results show that the sensitivity achievable by 
Planck is more than an order of magnitude better than 
that by current galaxy-clustering surveys and is compa- 
rable to that of larger next-generation surveys. 

This paper is organized as the following. In Sec. |TT] 
we first review the parametrization of the effect of fossil 
fields on scalar perturbations. Since this paper is con- 
cerned with the calculation of observables on a spherical 
sky, we recast the parametrization in Ref. Il6ll in terms of 
total-angular-momentum (TAM) waves jl8| . The TAM 
formalism makes the rotational symmetry of the observed 
sky manifest throughout and thus greatly simplifies our 
calculation. We then move on to present the calculation 
of the BiPoSHs and the BiPoSH power spectra. After 
that, we construct the minimal- variance quadratic esti- 
mator for the reduced amplitude, and calculate the cu- 
mulative signal-to-noise. We proceed to present numer- 
ical results for a local scalar-scalar-fossil bispectrum in 
Sec. mil They are complemented by Sec. IIVI which dis- 
cusses a real-space picture of the fossil-field effects on the 
temperature map. Finally, we make concluding remarks 
in Sec. El We also comment there on a divergence in the 
predicted observables from the scalar-scalar-tensor cor- 
relation in single-field slow-roll inflation. 



II. CMB BIPOSHS FROM INFLATION FOSSILS 

A. Fossil parameterization 

Interactions between a fossil field h and the curva- 
ture perturbation TZ can generate three-point functions 
of the type {TZTZh) at horizon crossing that then convert 
into a primordial bispectrum involving two scalar-metric- 
perturbation $ modes and a fossil-field mode after the 
end of inflation. A realization of the h field then locally 
induces departures from statistical homogeneity in the 
scalar autocorrelation function, which then appears as a 
correlation between different Fourier modes of the scalar 
perturbation. A general parameterization, 

($(ki)<i>(k2));,P(K) H^TrfSoik, + k2 + K)fP{h,k2,K) 

x<,(K)fc?fc^[/zP(K)]*. (1) 



is written down in Ref. [16[ to account for the possi- 
bility that the fossil field h can be scalar, vector or 
tensor field. Here ki and k2 are wavevectors of the 
two scalar-perturbation modes, K is the wavevector of 
the fossil-field mode, and hP(K) is the Fourier ampli- 
tude of the fossil field. The Dirac delta function reflects 
the translational invariance of the underlying physics. 
We introduce a generic bispectrum shape function 
fJl{ki,k2, K), which is related to the primordial scalar- 
scalar-fossil bispectrum through i3$$?t(ki, k2, K) = 
PaK)f^ik,,k2,K)el,iK)kfkl where P^K) is the 
power spectrum for the fossil field. Moreover, the sym- 
metric three- by-three polarization tensor e^j|(K) with 
p = 0, L,x,y, X geometrically distinguishes modula- 
tions of scalar, vector, and tensor type through the az- 
imuthal dependence when the configuration is rotated 
about the direction of K. The trace tensor e^^(K) cx 5ab 
and the longitudinal e^j(K) cx {KaKb-K^Sab/3)K-^ de- 
scribe a scalar fossil field. Two transverse- vectorial ten- 
sors e^'il'{K) cx Ki^a.w^^ satisfying K'^w^'^ = describe 
a transverse-vector fossil field. Likewise, two transverse- 
tensorial tensors e'j^^^CK) satisfying A''^e^j;^(K) = de- 
scribe a transverse-tensor fossil field. 



B. The TAM formalism 

Since angular observables defined on the two- 
dimensional sky are involved, we adopt the total-angular- 
momentum (TAM) wave formalism, recently developed 
in Ref. [l^l , to take full advantage of the rotational sym- 
metry of the problem from the very beginning. In the 
TAM formalism, scalar, vector, or tensor fields in three- 
dimensional space are expanded in terms of a complete 
set of spherical waves, as opposed to the conventional ex- 
pansion in terms of plane waves. These spherical waves 
are eigenfunctions of the Helmholtz equation with wave 
number K and are eigenfunctions of total angular mo- 
mentum and its third component with quantum num- 
bers J and M. Although the scalar case has been long 
known as the Fourier-Bessel expansion, the TAM formal- 
ism greatly simplifies calculations when vector or ten- 
sor fields have to be projected onto the two-dimensional 
sphere. In our case here, since three-by-three symmet- 
ric polarization tensors are used in Eq. ([1]), scalar, vector 
and tensor fossil fields can all be incorporated into a sym- 
metric traceless tensor field /lab(x). As the tensor field is 
expanded in terms of TAM waves, the longitudinal mode 
hjj^j (K) describes a scalar fossil field, the two divergence- 
free vectorial modes h^^j{K) and hj^j{K) describe a 
transverse-vector fossil, and the two divergence-free ten- 
sorial modes h'jfj (K) and h'j^j (K) describe a transverse- 
tensorial fossil. 

In terms of TAM coefficients ^im{k) and h'jj^j{K) for 
the scalar perturbation and the fossil field respectively, 



3 



the modulation analogous to Eq. ([T]) reads 



X (47r)3(-i) 



h+h + J 



kik2 



a,T,fci 



(iimi) 



(x) 



(2) 



for a = L, VE, VB, TE, TB respectively. Here Jx) 



and ^'"}A/')af,(x) are TAM wave functions for scalar and 
tensor fields. The overlaps of three TAM wave functions 
have been worked out in Ref. uM- Note that statistical 



homogeneity and isotropy require that fY^{ki, k2,K) 
flB{k^M,K) = fYAkiM.K) and fl^ik,,k2,K) 
frikuk2,K) = ff^{kuk2,K). 



C. CMB BiPoSHs 

Through the epochs of radiation domination and mat- 
ter domination to recombination, scalar metric pertur- 
bations source temperature and E'-mode polarization 
anisotropics of the CMB. Therefore, cross-correlations 
of different scalar-perturbation modes, as a consequence 
of modulation by the fossil field, give rise to cross- 
correlations between different harmonic modes of the 
CMB anisotropics. For clarity, in this paper we only 
discuss the effect on the temperature map. Still, the 
inclusion of E'-mode polarization is straightforward and 
will improve the overall sensitivity of detection, once the 
cross-correlations between temperature and E-mode po- 
larization are properly taken care of. 

For a given fossil-field configuration, the modulation 
violates statistical isotropy in the anisotropics. This ef- 
fect can be conveniently quantified in terms of a bipolar 
spherical harmonic (BiPoSH) expansion Being the 
two-sphere analog of the fossil-field modulation, Eq. ([Ij, 
it provides a parameterization for the most general two- 
point correlation function on the sky. Specifically, cross- 
correlations of CMB temperature multipoles read 



^(-f)"^(/iTOiZ2,-m2|JA/)Af^ 



JM 
i2' 



JM 



(3) 



where {Iimil2m2\l3m^) denotes the Clebsch-Gordan co- 
efficient. The BiPoSH coefficients ^fj^^, given by 



■7712 , — -A/ 



{aJ,,nA2*n2)h^ (4) 



exist if J, li and I2 can form a triangle. The angle brack- 
ets with a subscript h denote an average over all real- 
izations of $ for a fixed realization of h. Here we use 



V\^m^m2'm3 ^ compact notation for the usual Wigner-3j 
symbol. Since the fossil is parameterized by a symmetric 
tensor field, it only induces BiPoSHs with J ^ 2. 

Under the TAM basis, the scalar perturbation can 
be described by TAM coefficients <i>;„i(fc), which source 
CMB temperature multipoles according to 



f 

2^ 



(5) 



where gf{k) is the scalar radiation transfer function for 
temperature. 

Combining Eqs. g]), and (P, the CMB BiPoSH 
due to modulation of a single TAM wave of the fossil 
field is. 



aJM\ 



/(2/i + l)(2/2 + l 
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\ 47r 
y. f^{k^M.K)I^,i,j{kiM.K), 



kidkigi^{ki) I k2dk2gi^{k 



(6) 



where the parity P{a) = for a = L,VE,TE 
and P{a) = 1 for a VB,TB. The functions 
^hhj(ki,k2,K) are given by 



ir,i.Aki,k2,K) 



47r 



(2Zi + l)(2;2 + l)(2J+l) 



1/2 



(7) 



where the relevant overlaps of three TAM wave-functions 
can be found in Eqs. (80), (81), (84), (85), and (87) of 
Ref. [1^, for each a respectively. Note that the depen- 
dence on azimuthal quantum numbers toi, TO2, and M 
cancels out on the right hand side, and the final result for 
I;";^ j(fci, ^2, K) is reduced to an integral over the radial 
coordinate. 

In the above result for the BiPoSH coefficients, rota- 
tional invariance is manifest in that a given ^j'^^^ only 
generated by TAM waves of the fossil field with the same 
total-angular-momentum quantum numbers J and M. 
Moreover, due to parity conservation, L, VE, and TE 
modes only induce even-parity BiPoSHs, i.e. J + I1+I2 = 
even, while VB and TB modes only induce odd-parity 
BiPoSHs, i.e. J + li + I2 = odd. Therefore, vector 
and tensor fossils, both containing i?-mode TAM waves, 
can be distinguished from scalar fossils from their sig- 
nature in odd-parity BiPoSHs. Nevertheless, vector and 
tensor fossils cannot be geometrically distinguished from 
each other from CMB BiPoSHs, as they can with three- 
dimensional surveys [11] . This is heuristically understood 
as information loss when the three-dimensional correla- 
tion function of the scalar perturbation is projected onto 
the sky to give two-dimensional angular correlation func- 
tions of the CMB. 
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In addition to primordial mechanisms, late-time effects 
can distort a Gaussian, statistically isotropic map as well. 
Particularly, weak-lensing of the CMB by the foreground 
matter distribution also produces even-parity BiPoSHs, 
mimicking the effect of inflation fossils. Still, we antici- 
pate the spectra of BiPoSHs from lensing to differ from 
those from modulation by fossils. The reason is that the 
power spectrum for lensing by the scalar potential peaks 
at J ~ 60 [2l[ , but, e.g., for a local scalar-scalar-fossil bis- 
pectrum the fossil modulation effects dominate at J < 5, 
as our results will show later. This implies that the shape 
of BiPoSH power spectra can be used to break the degen- 
eracy. Besides, weak-lensing generates S-mode polariza- 
tion of the CMB, while a primordial scalar-scalar- fossil 
bispectrum does not. 



D. BiPoSH power spectra and estimators 

A stochastic background of the fossil field is presum- 
ably generated during inflation, just like the scalar field 
or inflationary gravitational waves, and it is characterized 
by a power spectrum, 

h}MiK)hj''M'{K'))-PHiK)^-^5niK-K') 

X SjJ'SMM'Saa' ■ (8) 

Statistical homogeneity and isotropy guarantee that 
Ph^'iK) = P^iK) ^ P^K) and P^^K) = 
Pr{K)^PUK). 

The statistical significance with which any individual 
BiPoSH is detected to be nonzero is expected to be small. 
In light of this, we average over all realizations of the 
fossil field, and define bipolar auto- / cross-power spectra. 



(jJ 




(9) 



to statistically measure the imprint of fossils. These 
correspond to four-point correlations in the temperature 
map. They are invariant under rotations. We empha- 
size that since the fossil-field power spectrum is statis- 
tically homogeneous and isotropic, statistical isotropy of 
the CMB is resumed after taking an ensemble average 
over the fossil field. 

To measure BiPoSHs from data, we use quadratic es- 
timators. Estimators for BiPoSH coefficients can be con- 
structed by 



J2 {~irHlim,h,m,\JM}al,^^al*„,^. (10) 



Then estimators for bipolar power spectra can be written 
down. 



1 



2 J 



T t 

M=-,J 



tl t2 



h+l2 + J\ 



(11) 



Note that the second term is needed to unbias the es- 
timators with respect to a Gaussian CMB map without 
fossil effects. 

We assume a phenomenological parameterization for 
the fossil-field background, which is described by two pa- 
rameters. One is the normalization for the power 
spectrum. 



p^{K)^r^p^{K), 



(12) 



with a fiducial-power spectrum shape P^{K). The other 
is the amplitude Bf^ of the scalar-scalar- fossil bispectrum. 



(13) 



with a fiducial bispectrum shape f^{ki,k2,K). Here Z 
can be scalar Z = L, transverse vector Z = {VE, VB} 
or transverse tensor Z = {TE, TB}. For vector or tensor 
fossils, both E modes and B modes have to exist. 

With these parameterization, Eq. ([6]) can be cast into 
the form 



A 



it\H^,,,(K)=^'hFi';ZiK)h%AK), (14) 



where a E Z and the coefficient function -f/^';" {K) can be 
read off from Eq. ([6|) as 

Pl^^iK) = -(-.)./(-l)'.+'^+n")i^ j^(2/, + l)(2/. + l) 

X J kldkigliki) J kidk2glik2)f^iki,k2,K) 
xir,i,j{k,M,K). (15) 
The BiPoSHs power spectra are then calculated to be 



r-' 



j^J: [ K^dKP^{K)F,';Z{K) [Fi'f^{K) 

(16) 



- ■^h-^hl2,hlV 



where we define the reduced amplitude A^ = ^^{B^)^ 
of the fossil background. There is an observational degen- 
eracy between the fossil-field power spectrum and scalar- 
scalar-fossil bispectrum, so only this combination is mea- 
surable. An estimator. 



A = IT 



(17) 



for the reduced amplitude A^ can be constructed, for 
each possible combination of J and i = 1, . . . , 4. Tak- 
ing into account the symmetry of ^^'^^^ j^^, without loss 
of generality we can fix h ^ h ^ h, and I2 ^ h if 

^1=^3- 

We can combine these estimators to bring down the 
statistical error. Under the null hypothesis, it can be 



shown that A 



's, being quartic expressions in af^, 
are nearly uncorrelated with each other. This is seen 



h ,lil2 ,131, 
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through calculating the null-hypothesis correlation ma- 
trix by applying the Wick expansion. Cross-correlation 
elements arc suppressed by either negative powers of 
{2li -I- 1) or a Wigner-6j symbol relative to the auto- 



We numerically evaluate and tabulate coefficient func- 



correlation elements. Treating Ai'f , , , 's as statisti- 
cally independent estimators, we then linearly combine 
them, and write down the inverse-variance-weighted es- 
timator. 



h,lil2.l3l4 




(18) 



where I2 h h) ^OOP^ over all independent combina- 
tions of multipoles. The denominator is equal to the 
inverse variance of , which is worked out to be 



Z\-2 _ 



2J + 1 



8 E/ E/ (jTT (jTT (jTT (jTT 
J hhhU '1 '2 (3 U 

J.Z IjrJ'Z ]* (19) 



X J", 



Hence, a high signal-to- noise S = A^ja^ indicates a 
detectable BiPoSH signature in the CMB from inflation 
fossils. 



III. NUMERICAL RESULTS 

Our discussion so far applies to the most general in- 
flation fossil with arbitrary power spectrum P^(K) and 
bispectrum shape f^{ki,k2, K). In order to provide some 
illustrative numerical estimates of the detectability of the 
signal, we specialize to the case where the bispectrum 
shape is of the local type. Such a bispectrum may arise 
if the fossil-curvature interaction is local. This occurs, 
for example, for the graviton-curvature-curvature bispec- 
trum 01 in single-field slow-roll inflation. We approx- 
imate the bispectrum in this case by the fiducial form 
f^{ki,k2,K) = (fcifc2)^'^/^ it attains in the squeezed 
limit fci,fc2 ^ K where it peaks. We assume a scale- 
free power spectrum (K) = 1/K^ for the fossil field 
as should arise if no physical scale other than the Hubble 
scale comes into play during inflation. 

Specializing to those fiducial forms for the power spec- 
trum and the bispectrum, we numerically compute sen- 
sivities a^^. We assume the standard flat ACDM cosmol- 
ogy with the WMAP+BAO-l-ifo best-flt cosmological pa- 
rameters taken from TABLE 1 of Ref . [22| , except that 
for consistency we assume perfect scale- invariance du ring 
inflation by setting = 1. The public code CAMB [23 
is used to tabulate radiation transfer functions. CMB 
multipoles up to /max = 3000 are considered. 



tions Fff^{K) from Eq. (HI]) for 5 x IQ-^Mpc"^ =^ iiT ^ 
5 X 10~^Mpc~^, ranging from fossil-field modes that are 
well outside of the horizon today to modes that were 
inside of the horizon at the last-scattering surface. Ex- 
amples of Pf^'i^ [K) are plotted in Fig. [2] to show some 
of the features. In the TAM-wave picture, the oscilla- 
tory nature of the F;'^;"(i^)'s can be understood as the 
value of the TAM wave function on the surface of last 
scattering. As K increases, spherical waves shrink in- 
ward to the origin, and hence nodes and anti-nodes al- 
ternate to pass the surface of last scattering. On the 
other hand, for large scale modes with very small values 
of K, even the first peak is well beyond the last scat- 
tering surface, and consequently BiPoSHs vanish in the 
infrared. Exceptions are the J = 2 even-parity modes 
a = L, VE, TE, whose wave functions do not vanish 
at the origin. Therefore, logarithmic infrared divergence 
arises in even-parity quadrupolar BiPoSHs due to super- 
horizon TAM modes. In this work, we introduce an in- 
frared cutoff at Kmin = 5 X 10~^Mpc~"'^, which does not 
affect the results except for the even-parity J = 2 Bi- 
PoSHs. 

We calculate sensitivities for Af^ at 3cr statistical signif- 
icance (Fig. [3]) as a function of the largest CMB multipole 
^max dominated by the signal. No instrumental noise is 
assumed for I ^ /max- We consider scalar, vector or ten- 
sor fossil fields Z = L,V,T respectively using two strate- 
gies. One is to fully exploit the information from both 
even- and odd-parity BiPoSHs. Since the signal-to-noise 
is dominated by even-parity BiPoSHs, this optimizes the 
sensitivity. In reality, however, lensing by gravitational 
potentials and other late-time mechanisms kick in on 
small angular scales / > 1000. They mimic the effect 
of fossil fields by generating even-parity BiPoSHs as well, 
making it necessary to disentangle between them. To cir- 
cumvent the problem, the second strategy is to include 
only odd-parity BiPoSHs. They provide clean probes of 
fossils with non-zero spin, with moderate compromise in 
the overall sensitivity. Even though there is degener- 
acy between the reduced amplitudes A]^ and Aj^ from 
measurement of the CMB BiPoSHs, one can still break 
it with a three-dimensional measurement from galaxy- 
clustering [l6j . 

The BiPoSH signal is dominated by anisotropic mod- 
ulation on the largest angular scales, and J ^ 5 almost 
saturates the detectability. On the other hand, modu- 
lation of small-scale anisotropics (very large li and ^2) 
contributes large signal-to-noise to bring down cr^ which 
is found to scale as ^maxj agreement with similar stud- 
ies of CMB trispectra [ij, [2^. Such a scaling law is 
justified analytically in the Sachs- Wolfe limit, valid on 
large scales. In that case, the radiation transfer functions 
are approximated by gf{k) = —ji{kr*)/3, where r* « 
14 Gpc is the comoving distance to the last scattering sur- 
face. For a given BiPoSH angular scale J and CMB scale 
li,l2 satisfying i ^ /i sa ^2 ~ I12 and L -^l^^ U~ I34, 
the signal squared can be shown to scale as ^^2^/34'^, while 
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Figure 2: Examples of coefficient functions F;^'j"(A') with A/ — h — h- Note that we plot —iFi^^{K) for J = odd since 

F;'^';"(isr) is imaginary in that case, (a) Top left panel: Pf^'ll of different polarization types a = L,VE,VB,TE,TB have 
different peak locations. Parity-even L, VE, TE modes contribute infrared divergence to the quadrupolar J = 2 BiPoSHs. 
Parity-odd VB,TB modes are not relevant to this infrared divergence, (b) Top right panel: as J increases, extrema and nodes 
of Ff^'^^ are systematically shifted to smaller scales, and the amplitudes decrease as well, (c) Bottom left panel: increasing the 

CMB multipole h rescales the amplitude of -Fj^'j", but extrema and nodes are marginally shifted, (d) Bottom right panel: the 
amplitude varies for different values of Al, and the locations of extrema and nodes change slightly. 



the noise squared damps out more rapidly as /]^2 '34 i so 



that each estimator A 



J,z 



contributes signal-to-noise 



squared oc ^12^34. Given the constraints \li — h] < J and 
1^3 — ^4! < J, the number of such estimators scales as 
■^Siia^i S;34^i • Therefore, the cumulative signal- 



E 



h2l 



12«34 OC 



to-noise squared is oc JJ2ii2<i„ 
^max- Nonetheless, at Zmax ^ 3000, gravitational lens- 
ing and instrumental error take over [12[ , cutting off the 
growth in sensitivity. 

To compare our results to the trispectrum parame- 
ter widely used in the local model of primordial 
non-Gaussianity [l3, [25l - l28| . we consider induced four- 
point correlations in (b) of Fig. [1] by integrating out 
an intermediate scalar fossil (dashed line), which mim- 
ics an intrinsic trispectrum. Since the induced trispec- 



trum is proportional to Aj^ , we find the conversion Aj^ = 
5.3 X 10~^^ ttvl- In this way, we compare our results to 
the WMAP 5-year bound \tnl\ < 3.3 x lO'' [H, which 
shows good discovery potential with Planck for both 
even-/odd-parity BiPoSHs combined and odd-parity Bi- 
PoSHs alone. We find good numerical agreement with 
the tnl forecasts of Ref. [ij]. The Planck satellite 
can improve upon WMAP's sensitivity to A^ by a fac- 
tor ^ 15. Planck's sensitivity corresponds to that of a 
galaxy-clustering surveys with volume V = 20 Gpc^ 
and a resolution /cmax = 0.1/iMpc~^. Thus, Planck's 
sensitivity is more than an order of magnitude better 
than that of the current SDSS-III BOSS survey [s^. 
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Figure 3: Predicted 3a sensitivity for the fossil-field re- 
duced amplitude Ah as a function of the maximum mul- 
tipole /max- Only BiPoSHs with J ^ 5 are considered, as 
they dominate the signal. For the quadrupolar signal J — 2 
from L, VE, TE modes, we cut off the infrared divergence at 
if = 5 X IQ-^Mpc-^ We relate our results to the primor- 
dial trispectrum parameter tnl in the local model of non- 
Gaussianity through A'h ~ 5.3 x 10~^^rjvL. The solid hori- 
zontal line is equivalent to the WMAP 5-year upper bound 
|t"jvz,| < 3.3 X 10*, and the dashed horizontal line corresponds 
to the projected Planck sensitivity at Zmax = 1500 from the 
CMB trispectrum. 



IV. LOCAL-DEPARTURE MODULATION IN 
REAL SPACE 

There exists, for the local-type scalar-scalar-fossil, a 
simple and illustrative real-space picture of the effect of 
the fossil field on primordial perturbations. To see this, 
we write, 

<i'(x) = (l + ^K,{^)V-'V'^V'^ $^(x), a e Z. 

(20) 

This reproduces the squeezed limit X ^ of Eq. ([T]) 
with the local-departure form Eq. (fT3|) . Here <I'g(x) is the 
gaussian field, and 7^$ is the dimcnsionless normalization 
of the scalar-perturbation power spectrum. Note that the 
relation between ^>(x) and $g(x) is local. 

For CMB multipoles I < 100, the temperature is de- 
termined primarily by the value of the potential $ at 
the surface of last scatter. Therefore, consider one TAM 
mode of the fossil field. The distortions to hot and cold 
spots induced by the fossil field are proportional to the 
value of the TAM wave function ^'(j5v/)ab(x) s^""" 
face of last scattering. This reduces to three-dimensional 
tensor field that lives on the surface of the two-sphere. 
This tensor field is a linear combination of the five tensor 
spherical harmonics Y^jm) L, VE, VB, TE, TB, 

as shown in Eq. (94) of Ref. [l3 where the comoving ra- 
dial distance r corresponds to the surface of last scatter- 



ing. We clarify here that although these tensor spherical 
harmonics are only functions of two-dimensional n, their 
tensor value at each n can have components purely per- 
pendicular (/? = TE, TB) or parallel (/? = L) to the line 
of sight, or even partially perpendicular (/3 = VE, VB). 




Figure 4: Local distortion of hot /cold spots as a result of mod- 
ulation from tensor harmonics. Intrinsic (dashed) and dis- 
torted (solid) isothermal contours are shown, (a) longitudinal 
harmonic y( jj\/)ai,(fi)i (b) vectorial harmonics y^j^/j^^ (n); 
(c) tensorial harmonics ^(jf/^f (n). 

Consider a patch in the sky, which is small relative to 
the typical variation scale of the fossil field, but large rel- 
ative to the scale of hot and cold spots. Then according 
to Eq. ((20l) the fossil field acts as a constant modula- 
tion factor across the patch. From Eq. (|20p . the depar- 
ture is proportional to the double gradient of $g(x), but 
only the components purely perpendicular to the line of 
sight can be seen. These then distort isothermal contours 
on the sky. This implies that local distortion patterns 
are distinct for each of the five types of tensor harmon- 
ics (Fig. S]): (1) The longitudinal harmonic y(jM)abi'^) 
is purely radial, so from a projected view onto the sky, 
isothermal contours expand or shrink in an isotropic way. 
As a result, only the contrast between hot and cold spots 
are enhanced or reduced, but the shape and the center re- 
main the same. (2) The vectorial harmonics (^) 
establish a preferred direction at given ii from a pro- 
jected view onto the sky. Isothermal contours arc shifted 
when they are perpendicular to that preferred direction, 
but are unchanged if parallel to it. The outcome is that 
spots are shifted parallel or anti-parallel to the preferred 
direction, and hot spots and cold spots shift in the same 
direction. (3) The purely transverse tensorial harmonics 
^(JM^^(^) define two principle directions perpendicu- 
lar to each other at given n from a projected view onto 
the sky. Isothermal contours always shift to the hot side 
along one principle direction, and shift to the cold side 
along the other principle direction. This leaves the center 
of the hot and cold spots unchanged, but induces ellip- 
ticity by elongating the spot along one of the principle 
directions and squashing it along the other. 

On the scale of the fossil-field modulation, the local 
distortion pattern varies from patch to patch (Fig. [5|). 
Given fossil-field angular-momentum quantum numbers 
J and AI, the full sky is divided into regions centered 
at zeroes of the harmonic. For modulation by the lon- 
gitudinal harmonic Y^j^j~^^f^{n), the hot/cold contrast is 
enhanced in one of the two adjacent regions and is re- 
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Figure 5: Modulation of isothermal contours by different ten- 
sor harmonics: (a) intrinsic contours without modulation; (b) 
longitudinal type; (c) vectorial i5-type; (d) vectorial B-type; 
(e) tensorial i5-type; (f) tensorial i3-type. 



duced in the other. For modulation by the vectorial 
E-mode harmonic Y(^M)ab^^)^ rings of hot/cold spots 
around the center of the region either expand away from 
it, or sink toward it, and an adjacent region has the oppo- 
site effect. For the B-mode harmonic Y^M)abi^)' rings of 
hot/cold spots systematically rotate around the center of 
the region, with opposite rotating directions in adjacent 
regions. Similarly, for modulation by the tensorial in- 
harmonic y(jff)ab(fi): hot/cold spots are elongated along 
the radial/tangential direction with respect to the center 
of the region. By contrast, for the i?-mode harmonic 
'^(jM)abi'^) the direction of elongation is rotated by 45°. 
Again, in adjacent regions the effect is the opposite. 

In the far-field limit, i.e. small fossil-field wavelength 
relative to the comoving distance to the surface of last 
scattering, the TAM wavefunction *(jM)ab asymp- 
totically proportional to the corresponding type of tensor 
harmonic Y(jM)ab(^)- ^^^^ case, longitudinal TAM 
waves mainly modulate hot/cold contrast, divergence- 
free vectorial TAM waves mainly induce displacement, 
and divergence-free tensorial TAM waves mainly gener- 
ate ellipticity. However, the interesting scales for a bis- 
pectrum of the local type are large scales, and a given 
TAM wavefunction typically projects to different types 
of tensor harmonics at last scattering and generates dif- 
ferent types of real-space distortions simultaneously. 



V. CONCLUSION 

In this paper, we have investigated the signature of 
inflation fossils in the anisotropics of the CMB temper- 
ature. To take into account the possibilities that the 
fossil can be a scalar, a vector or a tensor, we have used 
a generic parameterization for cross-correlations of two 
scalar perturbation modes as a consequence of the fossil. 
In harmonic space, non-zero BiPoSH coefficients are in- 
duced; in real space, the effects show up as modulations 



of temperature contrast and displacement and ellipticity 
of hot and cold spots from patch to patch across the sky. 
We have provided results for BiPoSH coefficients and Bi- 
PoSH power spectra for general inflation fossil scenarios. 

As an important example, we have numerically ex- 
amined the case of a local-type modulation, featur- 
ing a scalar-scalar-fossil bispcctrum dominated by the 
squeezed limit and a scale-free fossil power spectrum. By 
constructing the minimum-variance quadratic estimator, 
we have found that the combination = {B^)^ can 
be measured. Our numerical results show that modu- 
lations on large angular scales dominate the signal, and 
the sensitivity improves as ^~ax- For an experiment with 
Planck's sensitivity ^max 1500 and with a combined 
analysis of both even-parity and odd-parity BiPoSHs, 
we find detectability at A^/IO''^^ = 4.3, 2.3, 2.3 for 
Z — L,V,T respectively, at 3a significance. On the 
other hand, odd-parity BiPoSHs alone are clean probes 
of vector/tensor fossils, with sensitivities Af^ /10~^^ = 
12.8, 5.5 for Z = V,T respectively. These results are 
comparable to a sensitivity for the trispectrum parame- 
ter in the local model tnl ^ 2000. They are also signif- 
icantly better than the sensitivity estimated to be avail- 
able with current galaxy-clustering surveys. 

Even better signal-to-noise is expected if the i?-mode 
polarization is also considered (the scalar-scalar-fossil 
bispectrum wc discuss here does not generate or affect the 
S-mode polarization). In that case, the BiPoSHs arising 
from TT-correlation, i?i?-correlation and Tii^-correlation 
are correlated with each other. Therefore, a correct treat- 
ment must involve inverting the covariance matrix. We 
do not include the full analysis with polarization in this 
paper, leaving it to future work. Also, it will not be possi- 
ble to geometrically distinguish between scalar-, vector- 
and tensor- type fossil field from the CMB. A possible 
solution might be to scrutinize the shape of the bipo- 
lar power spetra. The shape is found to be degenerate 
between vector and tensor fossils for a primordial bispec- 
trum of local type, but it might be more discriminative 
for models with primordial bispectra of other shapes. In 
the latter case, fossil modes smaller than the scale of the 
last scattering surface can be important, and therefore 
distince BiPoSH shapes can be induced due to different 
polarization patterns for the scalar-, vector- and tensor- 
type fossil TAM wavefunctions in the far-field limit, as 
discussed in Sec. IIVI 

Weak lensing of the CMB by the foreground mass dis- 
tribution also generates non-zero BiPoSHs and could thus 
constitute a background for the fossil-field signal. How- 
ever, the fossil-field signal peaks (at least for a local- 
type coupling) at J < few, while the weak-lensing Bi- 
PoSHs peak at J ~ 100 and will have a very precisely- 
predictable shape. Vector and tensor fossil fields will 
moreover be distinguishable by the odd-parity BiPoSHs, 
which are not induced by lensing by density perturba- 
tions. Finally, lensing introduces a CMB-polarization B 
mode, which can further distinguish it from the effects of 
a fossil field. 
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Before closing we comment on the nature of the in- 
frared divergence in the prediction for the amphtude of 
the quadrupolar (J = 2) BiPoSHs that arises if the fos- 
sil field has a scale-invariant spectrum and an inflaton- 
inflaton-fossil three-point function of the local type. The 
predicted amplitude for the J = 2 BiPoSH then depends 
logarithmically on the smallest wavenumber K for the 
fossil field, and thus, on the onset of inflation. A simi- 
lar divergence arises also in the fossil-field prediction of 
galaxy clustering [l^. The result seems to imply that 
distance scales well beyond the horizon are having signif- 
icant effect on observables within the horizon. A similar 
logarithmic divergence has been discussed in Refs. f9l.[To| 
where it was argued that the tensor-scalar-scalar three- 
point correlation in slow-roll infiation may give rise to an 



power quadrupole that could be probed observationally 
jsHlll]- Ref. [i^l has argued, however, that the diver- 
gence is a gauge artifact. We believe it is important to 
fully understand this effect, as it occurs from the scalar- 
scalar-tensor three-point function that arises in single- 
field slow- roll inflation [3] , even without the introduction 
of additional new field. 
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